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We study the long time behaviour of a nonlinear oscillator subject to a random multiplicative 
noise with a spectral density (or power-spectrum) that decays as a power law at high frequencies. 
When the dissipation is negligible, physical observables, such as the amplitude, the velocity and the 
energy of the oscillator grow as power-laws with time. We calculate the associated scaling exponents 
and we show that their values depend on the asymptotic behaviour of the external potential and 
on the high frequencies of the noise. Our results are generalized to include dissipative effects and 
additive noise. 
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' The interplay of randomness with nonlinearity gives rise to a variety of phenomena both of practical and theoretical 
I significance that have been the subject of many studies since the pioneering work of Stratonovich [U, |3j H, 0] For 
■ example, electronic and Josephson junctions subject to thermal noise [^,0], the Faraday instability of surface waves 
and even human rythmic movements [8,] can be modelized by nonlinear oscillators subject to various sources of 
. noise (for a general reference on this subject, see e.g., @). In the regime where the dissipation of the effective nonlinear 

O ■ oscillator can be neglected, the random perturbation injects energy into the system and physical observables such as 
the oscillator's mechanical energy, root-mean-square position and velocity, grow algebraically with time fldl. fllL fT3|. 
ps| ■ This instability saturates when dissipative effects (that scale typically as the square-root of the energy) become large 
. enough to overcome the injection rate and the system reaches an non-equilibrium steady state. In series of recent 

' papers we have derived analytical results for quasi-Hamiltonian nonlinear random oscillators subject to 

external or internal random perturbation that is either a Gaussian white noise or an Ornstein-Uhlenbeck process. In 
particular, we have calculated the scaling exponents characterize the algebraic growth of the energy, the position and 
the velocity of the oscillator. 

(~| I In the present work, we consider the general case a nonlinear oscillator subject to a multiplicative noise with an 
• ^ ' arbitrary correlation function. We show that the physical observables display a scaling behaviour in the long time 
Ch . limit and we calculate the associated scaling exponents. The formulae obtained are valid for any type of noise whose 
power-spectrum decays as a power-law at large frequencies; in particular, the results previously derived for white and 
I Ornstein-Uhlenbeck noises appear as special cases of the general expressions that are derived in the present work. 
^ Finally, we briefly explain how our results can be generalized to include dissipative effects, additive noise, and how 
the scaling behaviour is modified when the noise power-spectrum decays decays exponentially at large frequencies. 
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'nI" ; I. THE DYNAMICAL EQUATIONS IN ACTION-ANGLE VARIABLES 

d 

We consider a nondissipative oscillator of amplitude x(t), trapped in a nonlinear confining potential U{x) and 
subject to a multiplicative noise S,{t) : 
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Our aim is to study the effect of the statistical properties of ^{t) on the long time properties of the dynamical variable 

We restrict our analysis to the case where the dominant term in the potential U when |a::| cxd is an even power of 
X. A suitable rescaling of x allows us to write 

U with n>2. (2) 

As the amplitude x{t) of the oscillator grows with time, the leading behavior oilA{x) for |a;| oo only is relevant and 
Eq.([T]) reduces to 

^x(t)+:.(t)2"-i=x(t)e(i). (3) 



We shall analyse the motion of the nonlinear stochastic oscillator following the method explained in J^. Defining 
the energy and the angle variables, 



E^h^ + l.^", and 0=7:^r"""^^. (4) 



2n ' ^ (2n)i/2n 7p 
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we transform the coordinates in phase space from position and velocity to energy and angle: 

X = ^^s„(0), (5) 

X = (2n)^i?^<(0) , (6) 

where the hypereUiptic function s„ is defined by the relation Sn{4>) — x/E^ . Using the auxiliary variable f2, that 
satisfies 

VL = {2n)'^ E"^ , (7) 
equation ^ is rewritten [l^, as a system of two coupled stochastic differential equations 

tl = (n-l)s„(0)4(0)f(t), (8) 

' (9) 



This system is rigorously equivalent to the original problem (Eq. ^ and has been derived without any hypothesis on 
the random perturbation term ^{t). It appears clearly in this formulation that the angle t/i is a fast variable whereas 
the action f2 is a slow variable : an effective stochastic equation for Q can thus be derived by adiabatic averaging over 
the angular variable. 

We now specify the statistical properties of the random perturbation ^{t). We shall consider a stationary Gaussian 
noise of zero mean value. A Gaussian process is entirely characterized by its auto-correlation function, defined as 

sit'^t) = {mm)- (10) 

In Fourier space, the power-spectrum of the noise is given by 

/ + 00 P + OD 

dtexp{iujt)S{t) ^ dtexp{iLut){^{t)C{0)) . (11) 
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For instance, if C(i) is a white noise of amplitude V, we have S{llj) — V. When ^{t) is an Ornstein-Uhlenbeck process 
of amplitude V and autocorrelation time t, we have 

Sit'-t) = ^e-\'-''y- and = — ^ . (12) 

2 T 1 + UJ^T^ 

In the white noise case, we derived in [l3l | the following scaling relations, valid when t — s- oo 

E ^ {Vt)^ , x--{Vt)^(^, x^{Vt)^^. (13) 

In the case when ^{t) is an Ornstein-Uhlenbeck process, we found in [l^ that, in the long time limit, the energy, the 
amplitude and the velocity of the oscillator grow as : 

E.(^)^. x^f^)^, .-(5)^^ (M, 



These colored noise scalings are always observed when t —^ oo, even if the correlation time r is arbitrarily small. 
We remark in equations (fT5|) and that the scaling exponents only depend on the asymptotic growth rate of the 
potential. However, the exponents are reduced by a factor 2 for Ornstein-Uhlenbeck noise as compared to their value 
for white noise; this means physically that the energy transfer from the noise to the oscillator is less efficient when 
the noise is correlated in time. In |14l |. we derived an analytic expression for the probability distribution function 
(PDF) in the long time limit that allowed us to derive not only the scaling laws but also to calculate all the prefactors. 
However, we did not find any satisfactory explanation that would explain the origin of the factor 2 reduction in of 
the scaling exponents. We show in the next section by analyzing the general case that the scaling exponents not only 
depend on the potential at infinity but also on the asymptotic behaviour of the spectral density of the noise at high 
frequencies. 



II. SCALING BEHAVIOUR IN PRESENCE OF ARBITRARY NOISE 



We consider now the more general case of a Gaussian noise which can be generated from the white noise by solving 
a linear differential equation of order a. Such a noise has a power-spectrum that decays as follows at high frequencies : 

- X>(wr)"2'" when |w|^oo. (15) 

The amplitude T) of the noise and the correlation-time r are defined by dimensional analogy with equation (|12p . The 
exponent a characterizes the high frequency behaviour of the the power-spectrum of the noise. We shall show in this 
section that, in the long time limit, the scaling behaviour of the nonlinear oscillator in presence of a noise with a 
power-spectrum that satisfies equation is given by 
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Dl ^^ 2(<T+i)(,i-i) 
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For (7 = (respectively a = 1) we recover the scalings given in equation (|13p for a white noise (respectively the scalings 
given in equation (|14p for Ornstcin-Uhlcnbeck noise). The formulae (|16p provide us with the explicit dependence of 
the growth exponents on the confining potential and on the noise spectral density. In particular, we observe that 
when the relative weight of the large frequencies is decreased (by increasing the value of a) the effective diffusion in 
the oscillator phase space becomes slower. 

We now derive the power laws p6p . One method to obtain these scalings is to apply the recursive adiabatic averaging 
procedure developped in [3, fl5l | . This advantage of this technique is to provide explicit expressions for the prefactors 
in the scaling law. The basic idea is the following : the naive averaging procedure applied to the system (|8l9p fails 
if the correlations between the angle </> and the noise ^{t) are discarded. However, by taking the derivatives of the 
equations ()8l9p repeatedly, a white noise contribution will appear (after a derivatives) and will provide the leading 
term in long-time behaviour of the oscillator. However, in order to apply this procedure, we need an explicit differential 
equation that relates the noise ^{t) to the white noise. Furthermore, this technique, though quantitatively precise, 
leads to very intricate calculations and is suitable for simple noises such as the Ornstein-Uhlenbeck process or the 
harmonic noise. 

We shall therefore derive here the power laws (fT6|) from a simpler argument. First, because we only want to extract 
the leading scaling behaviour when t — > 00, we simplify the equations ([5]) and ([5]), as we did in M : we keep only 
the dominant terms, replace the hyperelliptic function Sn{4>) by the circular function sin0 and put all the numerical 
factors to 1. We thus obtain the following system 

f2=C(i)sin0 and = , (17) 

which describes a pendulum with random frequency [17]. Now, the second order cumulant expansion [I61] of the 
stochastic Liouville equation for the PDF Pt{ft,<j>) associated with equations ^T7\ is given by 



-^ = LoPt + rf0(Li(t)exp(Lo0)Li(t-0)exp(-Lo(?))Ft 



(18) 



where the differential operators are defined as 



UPt = -g^i^Pt), (19) 

U{t)Pt = -A(^(i)sin0P,). (20) 
Following the method developed in [lB|, we evaluate the right hand side of equation (fT8|l by using to the formula 

exp(A)B exp(- A) = S + [yl, + i [A, [A, B]] + 1 [A, [A, [A,B]]] + ..., (21) 



with A — LiQ and B = Lii(t — 6). Wc find by induction tlie expression of the 71-th commutator 

[Lo, [. . . , [Lo, [Lo, Li(t - 0)]] ...]^at- 0) <t>) + ^^2"^(^' <t>) ] , (22) 

where the functions iJj""* and i?2"'' are given by 

H["\n,<l>) = (-l)"-if7"sin(0 + n|), (23) 
H^"\n,(l)) ^ cos{(f> + n^) . (24) 

Substituting the equations ()21ll24p in equation (|18p. we obtain, after taking the average over the noise and integrating 
over the variable 9 : 



J d9{U{t) exp(Lo^?)Li {t - 6) exp(-Lo^^)>P^ = sin f ^^i(^' *)^* + g^'^^^^^ ^)^* 



(25) 



with 



Hi{n,^,t) = y 1 \ ''^^ / d0r5((?)=- / sin(r!(? + 0)5((?)d0 , (26) 



n=0 

7^2(^^,0,^) = y^_L_i^ / d06'"5(6')= / sm{ne - (j))es(e)de . 
t'o Jo Jo 



(27) 



In the limit t 00, we average the differential operator (|25p over the fast variable and we deduce that 
r d6 f°° I d'^ f°° Id f°° 

where iS'(r2) represents the derivative of the power-spectrum with respect to the frequency. Finally, we apply this 
averaged operator to Pt(yi), the probability distribution function of the slow variable and we derive, in the long 
time limit, an effective averaged Fokker-Planck equation for Pt(0) : 



This equation only involves the power spectrum of the noise. If we substitute in this equation the large frequency 
behaviour (|15p of 5, we conclude by elementary dimensional analysis that : 

By using the definition ([71) of and the relations ([5]) and ^ between the amplitude x, the velocity x and the energy 
E, the scaling relations are derived. 

Our calculation can be generalized to the case of a weakly dissipative system, by including a linear friction term 
in equation ([T]) with dissipation rate 7. In this case, the oscillator's amplitude does not grow without bounds but 
saturates after a time of the order I/7. The system reaches a non-equilibrium steady state described by a stationary 
probability, which is in general not a Gibbs-Boltzmann distribution. In the case of a vanishingly small rate, the 
saturation time is large and the averaging technique can be still be applied. It can be shown, using the same method 
as above, that the stationary probability behaves as 

^stat(^)~exp(-C^^^ ^-^ ) whenii;-.cx3, (31) 



where C is a constant. The analytical expression of C and of the prefactors in -Pgtat '^^'^ ^^^^ be obtained by the 
method used here (a recursive adiabatic averaging similar to the one used for white and Ornstein-Uhlenbeck noise 
would be needed to go beyond scaling). 

We have asuumcd hitherto that the noise f {t) is generated from white noise by solving a linear differential equation 
of order a. Therefore, £^{t) is Gaussian and its power-spectrum S{lj) decreases as a power-law at larges frequencies. 
However, the derivation of the effective Fokkcr-Planck equation ([^H) does not rely on these assumptions and can still 
be performed when the correlation function of ^(t) has long-tails, for example, when S{oj) ^ exp(— t|cj|) (Lorentzian 
power-spectrum). We then obtain that 

ri^ilog-, (32) 

T T 

i.e., the amplitude of the system grows logarithmically with time. Similarly, if the power-spectrum of the noise is 
Gaussian, i.e., S{uj) ^ exp(— (tw)^) at high frequencies, the action variable fl^ - increases logarithmically with time. 
If the power-spectrum has a high-frequency cut-off at luq, if S(uj) = for lo > loq we conjecture fl saturates at 
large times. This fact is true, in particular, when the noise is a deterministic circular function for example in the case 
^(t) = cos(co'ot). 

Finally, we emphasize that the method we have used here can also be applied to a nonlinear oscillator subject to 
an additive noise, described by the following Langevin equation 

^xW = -^ + e(t), (33) 

where the power-spectrum of the noise ^(i) satisfies the equation Then, by calculations similar to those described 
above, the following scaling behaviour is derived (an alternative method that leads to the same results would be to 
use the Markovian approximation for the energy dynamics in the low friction limit, developed in [l9l. [20l|) : 



^-(t^) ^ --(t^J ^ --y-^) ■ (34) 

In particular, when the noise is white, we recover the fact that the energy grows linearly with time. 



III. CONCLUDING REMARKS 



A particle trapped in a nonlinear confining potential and subject to a random noise undergoes anomalous diffusion 
in phase space. In this work, we have found analytical expressions for the associated diffusion exponents, by using 
an action-angle representation of the equations of motion that allowed us to separate the slow variable from the fast 
variable and to derive an effective Langevin dynamics for the slow variable. Our results are derived for an arbitrary 
correlation function of the noise and are valid for multiplicative noise as well as for additive noise. We have shown 
that the high frequency components in the noise spectrum play a crucial role : the fastest the power spectrum of 
the noise decays at high frequencies, the slower is the diffusion. This fact has the following physical interpretation : 
the energy transfer from the external driving to the oscillator is optimal when the driving frequency is of the order 
of the natural frequency of the oscillator (e.i?., parametric resonance is maximal for driving at twice the natural 
frequency). However, because of the nonlinearity, the period T of the underlying deterministic oscillator decreases 
with the amplitude. Therefore, at large amplitudes, the energy transfer between the noise and the system involves 
higher and higher frequencies. If these frequencies have a small weight in the noise power spectrum, the amplification 
of the oscillator is less efficient. Our calculations put this intuitive reasoning on a quantitative basis and can be 
generalized to include dissipativc effects. In presence of dissipation, the system reaches a nonequilibrium steady state 
in which the stationary PDF is not the canonical Gibbs-Boltzmann distribution but a streched exponential. 
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